Stress Tensor and Bulk Viscosity in Relativistic Nuclear Collisions 
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We discuss the influence of different initial conditions for the stress tensor and the effect of bulk 
viscosity on the expansion and cooling of the fireball created in relativistic heavy-ion collisions. In 
particular, we explore the evolution of longitudinal and transverse components of the pressure and 
the extent of dissipative entropy production in the one-dimensional, boost-invariant hydrodynamic 
model. We find that a bulk viscosity consistent with recent estimates from lattice QCD further slows 
the equilibration of the system, however it does not significantly increase the entropy produced. 



INTRODUCTION 

The formation of strongly interacting matter at 
supranuclear energy densities has recently been studied 
in great detail in nuclear reactions at the Relativistic 
Heavy Ion Collider (RHIC). The analysis of the data col- 
lected in these experiments [l|, [2|, y, |j| has shown that the 
matter flows very rapidly at the moment of its breakup 
into free-streaming hadrons. The collective flow also ex- 
hibits a large anisotropy in non-central collisions, char- 
acterized by the "elliptic" flow parameter v-i. These ob- 
servations are commonly understood to imply that (i) 
the quark and gluon matter undergoes rapid equilibra- 
tion with thcrmalization times smaller than 1 fm/c [5j, 
and (ii) the quark gluon plasma (QGP) maintains low 
shear viscosity r\ not much larger than the conjectured 
lower bound 47T77kss = 8 where s is the entropy den- 
sity [fj. This has led to the claim that the quark gluon 
plasma formed at RHIC is the most perfect liquid known 
in nature [7J. 

The emerging picture still has some uncertainties. The 
predictions of hydrodynamic simulations, especially for 
the flow anisotropy in off-central collisions, depend on 
the assumed transverse density profile [9(. Moreover, 
the possible role of contributions to transverse flow from 
the pre-equilibrium phase of the reaction is not settled. 
Most hydrodynamic calculations assume that no trans- 
verse flow is present at the time of initialization, usually 
chosen in the range r = 0.5 — 1.0 fm/c after the onset 
of the reaction, although there are good reasons to be- 
lieve that transverse pressure gradients existing at ear- 
lier times will contribute to the generation of collective 
flow even if the parton momentum distribution is still 
anisotropic [ldj . 

Recently, it was pointed out that further complications 
could come from the bulk viscosity ( of quark and gluon 
matter near the QCD phase transition 



perature and weak coupling where £ ~ a 2 T 3 / In a J 1 [13( | 
while r\ ~ T z /(a 2 s In a J 1 ) [14|- The difference can be 
traced back to the near conformal invariance of QCD at 
high temperature. Near the pseudo-critical temperature 
T c , however, QCD is far from being conformal, as can be 
inferred from the large peak of the "interaction measure" 
(e - 3P)/T 4 at T c found in lattice QCD simulations [HJ]. 
Indeed, the ratio of the bulk viscosity over the entropy 
density, C/s, on the lattice (in the quenched approxima- 
tion) was recently found [161 ] to exhibit a narrow peak 
around T c of order unity with £/s 3> rj/s. The bulk vis- 
cosity can also be related directly to the trace anomaly 
of the energy momentum tensor, and an estimate can 
thereby be obtained from lattice results for the interac- 
tion measure [12]. 

The large spike of the bulk viscosity near T c immedi- 
ately raises several questions [12| which we want to ad- 
dress here. However, we also want to take a look at 
the larger picture. Ideal hydrodynamics requires com- 
plete thermalization of the matter. Isotropization of the 
pressure P and consistency with the equation of state 
P(e) are necessary conditions for equilibrium. Viscous 
hydrodynamics permits certain kinds of small deviations 
from equilibrium. In the local rest frame, shear stress 
7r y describes deviations from isotropy of the stress ten- 
sor while bulk stress H5 lJ measures deviations from the 
equilibrium equation of state. In the second-order formu- 
lation of viscous hydrodynamics 171 lilt > the deviations 
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viscosity can be neglected compared to shear viscosity 
in many systems in nature. This was also shown to be 
true in quantum chromodynamics (QCD) at high tem- 



of the stress tensor from its equilibrium form are not 
prescribed by external strains, such as flow shear or flow 
divergence, but can be independently given and only re- 
lax to the externally forced values over time. The fireball 
produced in a relativistic nuclear collision starts out with 
highly anisotropic particle distributions and therefore a 
large pressure anisotropy, which may not be related to 
the imprinted flow field. It is thus useful to explore how 
the hydrodynamical evolution responds to different ini- 
tial deviations from local equilibrium. 

It is not always obvious how large the deviations can be 
before the hydrodynamical approximation fails. In prin- 



ciple, viscous corrections extend the reach of hydrody- 
namics in relativistic heavy ion collisions to earlier times, 
but the gradient expansion of the stress tensor may break 
down [19l |. Similarly, a large bulk viscosity £ around T c 
could mean that the matter is driven far from equilib- 
rium around the phase transition and may even develop 
dynamical instabilities [20|. It has also been claimed 12 1 
that bulk stress could contribute significantly to entropy 
production, with possibly profound consequences for the 
hadronization mechanism at T c . Furthermore, one ex- 
pects the viscous bulk pressure II to be negative for an 
expanding system (II = —(d^u^ to first order in gradi- 
ents), implying that the effect of bulk viscosity is to slow 
down the expansion of the system. This should be most 
relevant in the longitudinal direction where the pressure 
is also reduced by the viscous shear stress. Hence, one 
expects the system to spend more time around T c than 
predicted by ideal hydrodynamics. This leads to yet an- 
other interesting prospect: Perhaps the evolution of the 
system is more sensitive to the equation of state near 
T c and the order of the phase transition than one would 
expect from ideal hydrodynamics where the relation of 
the phase transition to observables has been found to be 
rather obscure 
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Here, want to explore these questions in a hydrody- 
namic model with a simple space-time structure, but 
realistic equation of state, and bulk and shear viscosi- 
ties motivated by lattice QCD. Our framework is a one- 
dimensional boost-invariant fireball with translational 
and rotational symmetries in the transverse plane as first 
described by Bjorken 22j. We discuss the time evolution 
of the components of the pressure using second order hy- 
drodynamics. Our treatment is comparable to that of 
Baym [23| and Heiselberg and Wang [3, |25| . who ex- 
plored deviations from, and the approach to, equilibrium 
in a boost-invariant expansion in the framework of the 
relaxation time approximation to the Boltzmann equa- 
tion. We will also explore the consequences of several 
different assumptions about the initial longitudinal and 
transverse pressure. 

Let us rephrase our catalog of questions in the con- 
text of our hydrodynamic model: (i) How far from the 
equilibrium pressure P is the longitudinal pressure P z 
throughout the evolution and in particular at the phase 
transition? (ii) How close to isotropy is the pressure ten- 
sor, i.e. how large is \P± — P z \/P? (iii) What is the 
amount of entropy Sn produced by the bulk viscosity 
compared to contributions from shear viscosity 5$ dur- 
ing the lifetime of the fireball? 

We are aware that the simplified treatment with 1- 
dimensional expansion has several limitations. This ap- 
proximation does not permit us to explore effects related 
to the transverse expansion and the more rapid cooling 
coming from it. We also cannot treat the flow anisotropy 
in the transverse plane seen in off-central collisions. Since 
boost-invariant hydrodynamics tends to overestimate the 



time required for the matter to cool below T c it makes 
our study conservative in the sense that the effects of the 
bulk viscosity are underestimated, compared with those 
expected in a full three-dimensional treatment. 



VISCOUS HYDRODYNAMICS 

In the center of the fireball in a nuclear collision the 
viscous stress-energy tensor in the local comoving frame 
has the form [26|, |27|, [HI : 
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with the transverse and longitudinal pressure 



Px 

P, : = P 



p + n + -$ 



(2) 



Here P denotes the (isotropic) pressure in thermal equi- 
librium, $ and n denote the non-equilibrium contribu- 
tions to the pressure coming from shear and bulk stress. 
In particular, the traceless shear tensor in that frame 
takes the form ir zj = diag($/2, $/2, — $) consistent with 
the symmetries in the transverse directions. We refer the 
reader to Refs. 
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28l [29| for more details. 
At early times, $ and n will be given by the initial 
conditions for the stress tensor established by the mecha- 
nisms of energy and momentum deposition in the nuclear 
collision. Reflecting the dilution effect of the expansion 
on the local longitudinal momentum distribution of par- 
tons, the pressure components are expected to satisfy the 
ordering 



P z = T zz < T x 



jiyy = p 



(3) 



One can argue that the largest physically meaningful 
value of <I> at a time r is $ = 4(P + n), which corre- 
sponds to P z = —Px- Such a negative value of the lon- 
gitudinal pressure arises, for example, when the matter 
is completely in the form of coherent longitudinal fields 
at very early times [10, ISfJ]- However, such a configura- 
tion is very far from equilibrium and the hydrodynamic 
approximation is surely invalid. Once decoherence of the 
field is reached at a time Td ec the components of the phys- 
ical pressure should be positive: P Zl Px > 0. We want to 
explore various scenarios for the time evolution of $ and 
n, assuming initial values satisfying $ < P + H at the 
start time of the hydrodynamic evolution To > Td cc - 

In first-order (Navier-Stokes) dissipative hydrodynam- 
ics the bulk and shear stress are related to gradients in 
the system via the bulk and shear viscosities through 



n 



-Cd'X , 7T^ = 2r/V^u 



(m„"> 



(4) 



where w M = (cosh 77,0, Osinhr;) in our case is the ex- 
pansion velocity, r\ is the space-time rapidity, V^ — 
Qfj, _ u ^{u v d v ), and (. . .) indicates a projection orthogo- 
nal to u M , symmetrization of indices and removal of the 
trace. 

We follow the spirit of the (Israel-Stewart) theory of 
second-order dissipative hydrodynamics [1 81 ] by assuming 
that the actual bulk and shear stress have the freedom 
to relax to their first-order values at rates governed by 
relaxation times rn and r n . The equations governing 
the longitudinal expansion of the medium in our case are 



then given by [29j, |3l|, |3J, 133J : 
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It was recently pointed out by Baier et al. [33j that the 
terms in the square bracket in ([6]) are required in a the- 
ory with conformal symmetry. Conformal symmetry is 
approximately realized in QCD at high temperatures. 
Since these terms have not been studied quantitatively 
we will examine their influence on the evolution below. 
The entropy density s obeys the equation [34J ]: 
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dr 
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(8) 



ts = dS/(dydA) is the entropy per unit rapidity y and 
transverse area A. ts is constant for ideal hydrodynam- 
ics. Bulk and shear stress relax toward their Navier- 
Stokes values. The late time behavior of these values for 
1-dimensional boost-invariant expansion is then directly 
given by Eq. Q as 
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The larger the shear viscosity, the more anisotropic the 
pressure remains at late times. 

Relaxation times and viscosities are related by coeffi- 
cients /3o and 02 which are determined by the underlying 
theory: 



m = CA) , Tn = 2t7/3 2 



(10) 



Kinetic theory of massless partons predicts a value, P2 — 
3/(4P) [3lLl32|. This leads to a relaxation time which is 
roughly given by 



-(kin) _ 



2ttT 



(11) 



In conformal hydrodynamics a different behavior is ob- 
tained by matching the asymptotic form of a boost- 
invariant, longitudinally expanding thermal medium in 



the N = 4 supersymmetric Yang-Mills theory to a hy- 
drodynamic evolution 33J, I35J 



r (SYM) _ 



2 -In 2 
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(12) 



Obviously the value of the relaxation time from kinetic 
theory is about twice as long as the latter. We will 
test both values below. In absence of further reliable 
predictions, we will always assume the same relaxation 
time for the bulk stress as a function of temperature, 
Tn(T) = T-^iT). The coefficient Ai in © was determined 
for supersymmetric Yang-Mills theory to be [33j 
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VISCOSITIES AND INITIAL CONDITIONS 

We have already specified the equations of motion 
and our choice for the parameters tv, rn in the pre- 
vious section. We further assume that the matter is 
characterized by a minimal shear viscosity, i.e. we set 
Tj = t]kss = s/(4tt) . This is not in contradiction to lattice 
QCD results [36|], which lie close to the KSS bound. The 
pressure anisotropies found in our calculation can there- 
fore be considered as a lower bound. For the equilibrium 
equation of state we use recent lattice QCD results for 
two light and one heavy quark flavors [15j. We have pa- 
rameterized the reduced equilibrium pressure P/T 4 and 
interaction measure (e — 2>P)/T A for N T — 6 lattices from 
this reference. The critical temperature is T c = 196 MeV. 

For the bulk viscosity we explore several options. 
Our starting point is a recent calculation by Meyer in 
quenched QCD [161] . This calculation is not directly com- 
patible with our equation of state from unquenched QCD. 
To deal with this problem, we have here chosen to param- 
eterize the dimensionless ratio (,/ s as a function of the 
dimensionless ratio u = (e — 3P)/(e + P). This choice 
corrects for the shift in the critical temperature, and it 
softens the very steeply peaked behavior of £/s found by 
Meyer in quenched QCD. We call the bulk viscosity re- 
sulting from this fit (0 • We note that several mechanisms 
can contribute to the bulk viscosity [37|. In QCD, £/s is 
not only a function of the interaction measure but also 
of the speed of sound, the quark condensates, and the re- 
laxation time scale of the compression mode 12|, |38| . We 
also note that the statistical and systematic uncertain- 
ties of the existing lattice results for quenched QCD are 
quite large. We account for these combined uncertainties 
by running the hydrodynamic evolution for several val- 
ues of C which have been obtained from £0 by multiplying 
with a scaling factor cq: C — c cCo: and by changing the 
width of the peak near T c by a scale factor 1/af. 

We show the results of our fits in Fig.[TJ The kinematic 
bulk viscosity fitted from the lattice results as a function 
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FIG. 1: (Color online) Kinematic bulk viscosity C/s and kine- 
matic shear viscosity rj/s as a function of temperature T. rj/s 
(dash-dotted line, black) is determined by the KSS bound 
rj/s — l/(47r), C/ s (solid blue) is derived by a fit to the re- 
sults reported in [lg], using the conformal measure lo as scal- 
ing variable. This fit is denoted as Cf = aj = 1. The blue 
dashed curve shows a modified parametrization for C/s with 
double the peak height (c^ = 2) and half the peak width 
(Of = 2). Insert: parameterization of the reduced interaction 
measure (e — 3P)/T 4 taken from [15|] as a function of T. 



of temperature (shown by the solid (blue) curve) is com- 
pared with the kinematic shear viscosity from the KSS 
relation -qjs = l/(Air) (dash-dotted (black) line). The 
dashed (blue) curve shows a modified parametrization 
for C/s with double the peak height and half the peak 
width (c^ = a^ = 2). The bulk viscosity exceeds the 
shear viscosity for temperatures below 220 MeV. The in- 
sert shows the interaction measure fitted from the results 
of Cheng et al 0. 

For all simulations we choose a starting time tq = 0.3 
fm/c which is smaller than the equilibration times esti- 
mated from ideal hydrodynamics and is compatible with 
the expected decoherence time [39|, |40| of the initial gluon 
field. The initial energy density is fixed to be e(ro) = 50 
GeV/fm 3 which corresponds to an initial equilibrium 
temperature of roughly 400 MeV. We discuss three dif- 
ferent initial conditions: 

(i) "equilibration" sets n(r ) = $(ro) = 0. 

(ii) "1 st order" uses the values given by first-order vis- 
cous hydrodynamics, II (r ) = -C(T )/t , $(t ) = 
4?7(Xo)/(3to) where To is the initial temperature at 
to- 

(iii) "anisotropic" uses II(to) = — C(^o)Ao as in (ii), 
but fixes $(to) = P(to)+H(tq) for vanishing initial 
longitudinal pressure. 

We also remind the reader that we will run the hydro- 
dynamic evolution both with the conformal terms in Eq. 




FIG. 2: (Color online) (a) Relative bulk and shear stress, 
— H/P (solid blue) and $/P (dashed black), as functions of 
time r for scenario (iii,S,SYM,c^ = 1), i.e. anisotropic ini- 
tial conditions, standard evolution without conformal terms 
and short relaxation times. The time when T c is reached is 
indicated by the triangle, (b) Relative transverse and longi- 
tudinal pressure, P±/P (red) and P z /P (black), as functions 
of time t for (iii,S,SYM,Cf = 1) (solid lines) and for the same 
scenario but with C£ set to zero (dashed lines), (c) Relative 
entropy production from bulk and shear stress, Sn/S/ (solid 
blue) and S<s>/Sf (dashed black) as functions of time r for the 
same scenarios with cq — 1. All results are for a^ — 1. 



^ (we will denote this scenario by the label "C" ) and in 
a standard version without them (denoted by "S"). We 
also explore both the short ("SYM") and long ("kin") 
relaxation time. Each run below will be denoted by a 4- 
component label indicating the set of initial conditions, 
relaxation times, absence or presence of the conformal 
terms and the scaling variable cq = C/Co- 



RESULTS 

In this section we present the results of numerical so- 
lutions of the viscous hydrodynamical equations (J3H7J)- 
We start by checking the influence of different values of 



the bulk viscosity by varying cq with a^ fixed. We com- 
pare vanishing bulk viscosity (cq = 0) with the lattice 
inspired value (cv = 1)- Fore these runs we fix the set of 
initial conditions to anisotropic (iii), use evolution with- 
out the conformal terms (S) , and choose the short relax- 
ation time (SYM). Figure [5] shows the development of 
different components of the pressure relative to the equi- 
librium pressure P as a function of time. Part (a) shows 
the dissipative stress components —H/P and <&/P due 
to bulk and shear viscosity, respectively, for Cf = 1. As 
expected the effects of shear viscosity are dominant at 
early times while bulk stress peaks later, when £/s starts 
to rise around T c . The peak in — H/P is reached already 
slightly before the critical temperature T c — the time 
t c w 8.4 fm/c when the system reaches T c (for cq = 1) is 
indicated by the solid triangles on the abscissae. 

Figure [2] (b) displays the relative longitudinal and 
transverse pressure P±/P and P z /P, respectively for 
cq = \ (solid lines) and cq = (dashed lines). Both 
quantities develop a pronounced minimum just above T c 
if bulk viscosity is present, which corresponds to the max- 
imum in — II/P. The minima above T c are absent for van- 
ishing bulk viscosity. We observe that the system is still 
very anisotropic even at T c , and with cq = \ both pres- 
sure components are well below the equilibrium value. 

As mentioned earlier, the stability of boost invariant 
relativistic hydrodynamics in the presence of a bulk vis- 
cosity was recently studied by Torrieri and Mishustin 
[20(. Their treatment differs from ours in two impor- 
tant details: They study the first-order (Navier-Stokes) 
formulation of viscous hydrodynamics, and they use a 
parametrization of the temperature dependence of the 
bulk viscosity which is much more strongly peaked near 
T c and attains a much higher peak value. Even for their 
most conservative choice (zq = 0.1 in the notation of 
[20(), the peak in £/s is about twenty times higher than 
our peak value and five times as high as the largest value 
obtained by Meyer [36( . We have solved our set of equa- 
tions with their parameters and used a very small relax- 
ation time Tn = r^ to emulate the Navier-Stokes limit. 
We find that P z becomes strongly negative in the range 
where (/s peaks, indicating that the matter is not only 
hydrodynamically unstable, as found in ref. [20J |. but also 
thermodynamically unstable. This behavior is obviously 
a result of the highly peaked parametrization adopted by 
the authors of ref. [20| • 

We now return to our own study. Figure [2] (c) shows 
the entropy per unit rapidity and transverse area pro- 
duced by shear and bulk viscous effects, S$ = rsq, and 
SV = rsn , respectively. They correspond to the first and 
second term in Eq. ([8]) and are shown relative to the total 
final value Sf = Tfs(rf) where the final time r/ is fixed 
at 50 fm/c. We note that despite its dramatic effects 
on the longitudinal pressure the contribution of the bulk 
stress to entropy production is rather moderate. The en- 
tropy produced by shear stress is much larger, due to the 
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FIG. 3: (Color online) (a) Relative longitudinal pressure 
P z /P as function of time r for scenario (i,S,SYM) and var- 
ious choices for the bulk viscosity as shown in the legend 
(c,; = 0,1,2 and a^ = 1,2. The solid line shows our "stan- 
dard" parametrization (cj = a^ — 1). (b) Relative entropy 
production from bulk stress, Sn/Sf as function of time r for 
the same scenarios as in part (a). 



large velocity gradient in the initial state which gener- 
ates large dissipative effects. However, the majority of 
the entropy production is confined to the earliest time 
period tq < r < 1 fm/c, suggesting that a hydrody- 
namical description of the matter rapidly loses reliability 
before 1 fm/c. We also remind the reader that our results 
are obtained for the minimal value of the shear viscosity 
(?/ = s/(4tt)). Had we chosen a larger shear viscosity, the 
produced entropy would be respectively larger. 

In Fig. [3] we have varied the parameters c^ and a^ , in- 
fluencing the height and width of the peak in the bulk 
viscosity near T c , respectively, for the scenario (i,S,SYM). 
The short-dashed (blue) curve in Fig. [3]Ja) shows the rel- 
ative longitudinal pressure in the absence of any bulk vis- 
cosity for comparison. The solid and dash-dotted (red) 
lines are obtained for our parametrization of C(^) as 
shown by the solid (blue) line in Fig. Q] with its height 
scaled by the factor c^ = 1,2. Obviously, the relative 
longitudinal pressure drops to almost zero above T c for 
cq = 2 indicating that the local equilibrium assumption 
begins to break down in this case. The dotted and long- 
dashed (black) curves correspond to the parametrization 
(a^ = 2) for a narrower peak in £(T), as shown by the 
dashed (blue) line in Fig. Q] In this case the onset of the 
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FIG. 4: (Color online) (a) Relative transverse and longitudi- 
nal pressure, P_l/P (dashed) and P z /P (solid), as functions of 
time r for (S,SYM,c^ = 1) and initial conditions (i) (black), 
(ii) (blue), and (iii) (red), (b) Relative entropy production 
from bulk and shear stress, Sn/S/ (solid) and S<s>/Sf (dashed) 
as functions of time r for the same set of scenarios. Note that 
the different curves for Sn/S/ lie almost on top of each other. 
The triangle indicates the time of critical temperature for ini- 
tial condition (iii). For conditions (i) and (ii) T c is reached 
slightly earlier. All curves are for a$ = 1. 



strong reduction in P z j P is delayed compared with the 
case olq = 1, but the effect in the immediate vicinity of T c 
(the location is indicated by the (green) triangle on the 
abscissa) is found to be mainly sensitive to the height of 
the peak, parametrized by c^, not its width. It is impor- 
tant to note, however, that the hydrodynamical evolution 
becomes increasingly sensitive to the precise value of the 
equilibration time for the bulk viscosity, rn, as the width 
of the peak in C,(T) becomes narrower. The small value 
chosen here (rn — r ff ) may be inappropriate for a very 
narrow peak. Insofar as a narrow, high peak in C(T) is 
indicative of a near-critical behavior of the medium near 
T c , one would expect any mode that participates in this 
behavior to exhibit critical slowing down and its relax- 
ation time to increase. 

We now compare the impact of different initial con- 
ditions on the standard evolution (S) for a bulk viscos- 
ity given by the scaling factors cq — a^ — 1 and short 
relaxation time (SYM). Figure [4] (a) shows the relative 
pressure components P z /P and P±/P. The most notice- 
able feature here is that even for equilibrium initial con- 
ditions (i) the strong gradients in longitudinal direction 
drive the system immediately off equilibrium. The first 



> 


0.5 


Q. 




O 


04 






H-» 




c 


3 


LU 




^^ 


0.2 


CD 




DC 


0.1 




0.0 




CD 



(b) 



Sn/S, 



___ S$/Sf # *"" 



1 - 4 ^:>. 

1.2 **::>^.\, 



— t„ kin. 

— t, conf. 

— t„ EOM conf. 

— t„ EOM conf., C=0 




FIG. 5: (Color online) (a) Relative transverse and longitu- 
dinal pressure, P_l/P (dashed) and P z /P (solid), (b) Rela- 
tive entropy production from bulk and shear stress, Sn/S/ 
(solid) and S<s>/Sf (dashed). The different scenarios shown 
for initial conditions (iii) are: r-n = m = t4 from kinetic 
theory, see Eq. {TT}, no conformal terms in <j6j , C( = 1 [red]; 
TV = tyi = tv from conformal symmetry, see Eq. (|12|l . 

no conformal terms in (|S}, c^ = 1 [blue]; tv — rn = tv , 

conformal terms in <(6j switched on, c^ = 1 [black]; r^ = 
tyi — t^ , conformal terms in (|6]) switched on, cq — 
[gray]. The triangles indicate the largest and smallest times 
of critical temperature corresponding to the last and the first 
scenarios, respectively. 



order initial conditions (ii) are very close to the maximal 
anisotropic initial conditions (iii) which were loosely ex- 
trapolated from classical gluon fields. Interestingly, the 
effect of the different initial conditions is wiped out after 
a very short time At « 0.5 fm/c and the system evolves 
in an universal way from that time forward. However, 
the different initial conditions for $ leave a trace in the 
entropy produced during this stage of the evolution, as 
can be seen in Fig. Q] (b) . On the other hand, the entropy 
production from bulk stress picks up most contributions 
around T c and is independent of the initial conditions. 

Lastly, we study the influence of the relaxation times 
and the conformal terms in the evolution. Figure [5] (a) 
shows P z /P and P±/P, using the standard equation of 
motion (S) for <1> and cq = a^ = 1 using both the estimate 
for relaxation times t w and Tn from kinetic theory (kin) , 
and the lower estimate (SYM). We also show the lower t„. 
(kin) with the additional conformal terms (C) switched 
on and c^ = 1, and the same with cq = 0. Figure O (b) 
shows the relative contributions to entropy production in 



these four scenarios as before. 

Obviously, larger relaxation times lead to increased en- 
tropy production and larger deviations from equilibrium. 
This is very clear for the shear contributions at early 
times. The effect of different relaxation times seems to 
be very much suppressed for the bulk stress. In fact, after 
r w 2 fm/c we again see a universal evolution of the sys- 
tem for fixed bulk viscosity. As expected the additional 
conformal terms in ^ lead to a suppression of the shear 
stress which manifests itself in a rapid relaxation away 
from the maximum anisotropic initial condition. Enforc- 
ing conformal symmetry leads to smaller anisotropies be- 
tween the transverse and longitudinal direction and the 
system is generally closer to equilibrium. This agrees 
with the recent observation of Song and Heinz 4l| made 
in the context of a study of two-dimensional boost invari- 
ant hydrodynamics including transverse expansion. 



DISCUSSION AND SUMMARY 

We can now answer some of the questions posed at 
the beginning of the paper. First, independent of the 
initial conditions we find that the longitudinal pressure 
reaches at most half of the equilibrium pressure through- 
out the entire lifetime of the quark-gluon plasma phase 
if the bulk viscosity is close to the values suggested by 
lattice QCD (cq rj 1), even if the shear viscosity takes its 
minimum value ??kss- This keeps the system away from 
equilibrium and reduces the work done in the longitudi- 
nal expansion. Bulk stress is the main contribution to 
this effect after about 2 fm/c and we find that the evo- 
lution of the system after this time, for a given equation 
of state at equilibrium, is solely determined by the value 
of the bulk viscosity and largely independent of the ini- 
tial conditions and relaxation times. Scaling our lattice 
inspired bulk viscosity with factors cq > 1 leads to un- 
acceptably small or even negative longitudinal pressure 
which would indicate a breakdown of the hydrodynamic 
picture. These results remain qualitatively unchanged if 
the width of the peak in the kinematic bulk viscosity is 
varied by a factor of two. 

We also found that isotropization of the stress tensor 
proceeds rather slowly with \Px — P z \/P a s large as 20% 
around T c . It was recently pointed out by Martinez and 
Strickland [42| that the anisotropy of the stress tensor of 
the quark-gluon plasma may be observable via changes 
in the dilepton yield. On the other hand, the strong re- 
duction in the relative longitudinal pressure, which may 
even lead to negative values of P z near T c , is reminis- 
cent of a first-order phase transition, where the negative 
pressure is avoided by the formation of a mixed phase. 
Under favorable conditions, the delay of the expansion 
caused by mixed phase formation can be observed as a 
directional dependence of the identical-particle correla- 
tion function in density interferometry 43|, 144 , I45J . The 



scenario found here may show similar effects, but differs 
from the traditional one by the anisotropy of the stress 
tensor caused by the continued presence of the shear vis- 
cosity. A realistic exploration of the influence of the bulk 
viscosity on identical particle correlations will the require 
the hydrodynamical treatment of the transverse expan- 
sion including both, bulk and shear viscosity. 

On the other hand, we find that bulk stress has a rather 
modest impact on entropy production, contrary to some 
previous expectations. For the parametrization C(T) con- 
sidered here (see Fig.[TJ c^ = 1) the entropy increase due 
to the bulk viscosity is at most 10% of the final entropy. 
The reason is that the velocity gradients due to the longi- 
tudinal expansion have decreased significantly when the 
fireball approaches T c , This can also be easily seen by 
noting that we have only shown bulk and shear stress rel- 
ative to the equilibrium pressure. While those are com- 
parable, by the time the maximum in —H/P around T c 
is reached the pressure P has dropped significantly and 
the absolute values of II are much smaller than the values 
of i> reached at times smaller than 1 fm/c. A bulk vis- 
cosity much larger than indicated by our extrapolation 
of the existing lattice QCD results (represented here by 
scale factors Cf 3> 1) would be necessary to dominate 
entropy production after decohercncc. This statement is 
independent of the initial condition for the bulk stress. 

Non- linear terms in the evolution of the shear stress, 
dictated by conformal symmetry, suppress the shear 
stress and lead to reduced anisotropies and entropy pro- 
duction. However, these effects are not large enough to 
affect the conclusions drawn above qualitatively. 

To summarize, large bulk viscosities around T c lead 
to prolonged deviations from equilibrium that could be 
sizable throughout the entire lifetime of the quark gluon 
plasma. Bulk viscosities just slightly larger than cur- 
rently favored could easily lead to a breakdown of the 
hydrodynamic approximation around T c . The decreased 
pressure should slow down the expansion of the system 
and increase the time spent in the vicinity of the phase 
transition. However, the amount of entropy produced 
through bulk stress around T c is smaller than that pro- 
duced by shear stress at earlier stages of the evolution 
and thus does not result in a large increase of the final 
particle multiplicity, unless the bulk viscosity is much 
larger than that considered here. 

Acknowledgments: We thank the members of the 
Yukawa Institute of Theoretical Physics for their hos- 
pitality during the YIPQS Molecule Entropy Production 
Before QGP. RJF would like to thank P. Huovinen for 
helpful discussions. This work was supported in part 
by the U. S. Department of Energy (grants DE-FG02- 
05ER41367, DE-AC02-98CH10886), the RIKEN-BNL 
Research Center, the Texas A&M College of Science, and 
the Bundesministerium fur Bildung und Forschung. 



[1] BRAHMS, I. Arsene et al, Nucl. Phys. A757, 1 (2005). 
[2] B. B. Back et al, Nucl. Phys. A757, 28 (2005). 
[3] STAR, J. Adams et al, Nucl. Phys. A757, 102 (2005). 
[4] PHENIX, K. Adcox et al, Nucl. Phys. A757, 184 (2005). 
[5] U. W. Heinz and P. F. Kolb, Nucl. Phys. A702, 269 

(2002). 
[6] P. Kovtun, D. T. Son, and A. O. Starinets, Phys. Rev. 

Lett. 94, 111601 (2005). 
[7] T. Hirano and M. Gyulassy, Nucl. Phys. A769, 71 

(2006). 
[8] M. Cheng et al, Phys. Rev. D74, 054507 (2006). 
[9] T. Hirano, U. W. Heinz, D. Kharzeev, R. Lacey and 
Y. Nara, Phys. Lett. B 636, 299 (2006). 
[10] R. J. Fries, J. Phys. G34, S851 (2007). 
[11] D. Kharzeev and K. Tuchin, arXiv:0705.4280 [hep-ph]. 
[12] F. Karsch, D. Kharzeev, and K. Tuchin, Phys. Lett. B 

663, 217 (2008). 
[13] P. Arnold, C. Dogan, and G. D. Moore, Phys. Rev. D74, 

085021 (2006). 
[14] P. Arnold, G. D. Moore, and L. G. Yaffe, JHEP 11, 001 

(2000). 
[15] M. Cheng et al, Phys. Rev. D 77, 014511 (2008). 
[16] H. B. Meyer, Phys. Rev. Lett. 100, 162001 (2008). 
[17] I. Miiller, Z. Phys. 198, 329 (1967). 
[18] W. Israel, Annals Phys. 100, 310 (1976); W. Israel and 

J. M. Stewart, Ann. Phys. 118, 341 (1979). 
[19] M. Lublinsky and E. Shuryak, Phys. Rev. C 76, 021901 

(2007). 
[20] G. Torrieri and I. Mishustin, arXiv:0805.0442 [hep-ph]. 
[21] D. H. Rischke and M. Gyulassy, Nucl. Phys. A608, 479 

(1996). 
[22] J. D. Bjorken, Phys. Rev. D27, 140 (1983). 
[23] G. Baym, Phys. Lett. B138, 18 (1984). 
[24] H. Heiselberg and X.-N. Wang, Phys. Rev. C53, 1892 



(1996). 
[25] H. Heiselberg and X.-N. Wang, Nucl. Phys. B462, 389 

(1996). 
[26] D. Teaney, Phys. Rev. C68, 034913 (2003). 
[27] A. Muronga and D. H. Rischke, (2004) arXiv:nucl- 

th/0407114. 
[28] A. Muronga, Phys. Rev. Lett. 88, 062302 (2002), 

[Erratum-ibid. 89, 159901 (2002)]. 
[29] U. W. Heinz, (2005), nucl-th/0512049. 
[30] Y. V. Kovchegov and A. Taliotis, Phys. Rev. C76, 

014905 (2007). 
[31] A. Muronga, Phys. Rev. C69, 034903 (2004). 
[32] R. Baier, P. Romatschke, and U. A. Wiedemann, Phys. 

Rev. C73, 064903 (2006). 
[33] R. Baier, P. Romatschke, D. T. Son, A. O. Starinets, and 

M. A. Stcphanov, JHEP 0804, 100 (2008). 
[34] U. W. Heinz, H. Song, and A. K. Chaudhuri, Phys. Rev. 

C73, 034904 (2006). 
[35] M. Natsuumc and T. Okamura, (2007) arXiv:0712.2917 

[hep-th]. 
[36] H. B. Meyer, Phys. Rev. D76, 101701 (2007). 
[37] K. Paech and S. Pratt, Phys. Rev. C 74, 014901 (2006). 
[38] S. Jeon, Phys. Rev. D 52, 3591 (1995). 
[39] B. Miiller and A. Schafer, Phys. Rev. C 73, 054905 

(2006). 
[40] R. J. Fries, B. Miiller and A. Schafer, arXiv:0807.1093 

[nucl-th] . 
[41] H. Song and U. W. Heinz, arXiv:0805.1756 [nucl-th], 

Phys. Rev. D (in press). 
[42] M. Martinez and M. Strickland, arXiv:0805.4552 [hep- 
ph]. 
[43] S. Pratt, Phys. Rev. D 33, 1314 (1986). 
[44] G. Bertsch, M. Gong and M. Tohyama, Phys. Rev. C 37, 

1896 (1988). 
[45] G. F. Bertsch, Nucl. Phys. A 498, 173C (1989). 



